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Abstract 

We define a (chaotic) deterministic variant of random multiplicative 
cascade models of turbulence. It preserves the hierarchical tree struc- 
ture, thanks to the addition of infinitesimal noise. The zero-noise limit 
can be handled by Perron-Frobenius theory, just as the zero-diffusivity 
limit for the fast dynamo problem. Random multiplicative models do 
not possess Kolmogorov 1941 (K41) scaling because of a large-deviations 
effect. Our numerical studies indicate that deterministic multiplicative 
models can be chaotic and still have exact K41 scaling. A mechanism is 
suggested for avoiding large deviations, which is present in maps with a 
neutrally unstable fixed point. 
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1 Introduction 



One popular way to describe the small-scale activity of fully developed turbu- 
lence is to suppose that energy is transferred from the injection scale to the 
viscous scales, by a multi-step process along the inertial range. This idea has 
been often used to predict important features of turbulent flows. Still, relations 
with the structure of Navier-Stokes equations are poorly understood. 

In his 1941 work, Kolmogorov [|I|] uses his phenomenology to postulate that 
all the statistical properties of the turbulent flow at scales belonging to the 
inertial range are universal, in the sense that they depend only on the scale £ 
and on the mean energy dissipation rate per unit mass e. As a consequence, 
moments of velocity increments, over small distances i, should posses univer- 
sal forms. This led Kolmogorov to dimensionally-based expressions for the 
structure functions : 

Si"^ = E{6v^,{x)} = E{{v{x + £)- v{x)y} ~ {e£y/\ (1) 

where E{- ■ ■} denotes ensemble average and the symbol ~ means equality 
within 0(1) multiplicative constants. Thus ([^) predicts scaling behavior for 
the structure functions, the function of order p having the exponent (p = p/3. 

Although early experimental data seemed to confirm the prediction (at least 
for p = 2), the K41 theory has been criticized because it does not take into 
account the natural (and also experimentally verified) presence of fluctuations 
in the energy dissipation. These fluctuations are commonly believed to be 
the consequence of the chaotic transfer of excitations along the inertial range. 
A possible way to account for this effect, suggested by Obukhov 0, is the 
following. First, one introduces the average rate of energy dissipation over a 
cubic box Af {x) of side i centered on x : 

= i , A^)^^- (2) 

Second, fluctuations in the velocity increment are related to fluctuations in 
ei{x) by a "bridging relation" suggested by Kolmogorov's 1941 theory, viz 

6vfix) ~ (q(x))V3. (3) 

This bridging relation is still widely used. It has received good experimental 
support P|. Still it leads to some conceptual difficulties which are unrelated to 
the material to be presented in this paper From (^ and (^, we find 

S^''^ ~ E{e//^} (4) 
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Considerable experimental @, §|, and theoretical work ^, ^ has been 
devoted in recent years to measuring and predicting the dependence of the 
quantities E{e^^^} on the scale i. To date, it seems that a certain consensus 
has been reached. From the experimental side, there is clear evidence of a 
non-trivial ^-dependency (at least for high order moments). Similarly, from a 
theoretical point of view, all models based on random multiplicative cascades 
(Section 2) introduce power-law corrections to K41 scaling. The fact that these 
simple models have deviations to K41 scaling has led perhaps to the miscon- 
ception that any cascade model having non-trivial fluctuations is inconsistent 
with K41 scaling. 

Actually, we shall show that the chaotic transfer of energy described by a 
deterministic multiplicative process can still be consistent with K41, or more 
precisely that possible deviations disappear in the "fully developed limit", i.e. 
for cascades with very many steps. 

The paper is organized as follows. In Section 2, we shall connect the well- 
known case of random independent multiplicative cascade models and our (de- 
terministic) case in two steps. First, we introduce a Markov random model and 
then, by means of a special limiting construction, we obtain the deterministic 
model. A particular class of deterministic multiplicative models is introduced 
in Section 3 and studied numerically in Section 4. Finally, in Section 5 we give 
a possible interpretation of the lack of corrections to K41 scaling observed in 
our model. 



Pure-random and noisy-deterministic multi- 
plicative models 



Random multiplicative models were introduced by Novikov and Stewart |13 



and by Yaglom as a simple way to describe stochastic transfer of energy 
along the inertial range. Their fractal properties were discussed by Mandelbrot 
rH|. Let us give now the definition of these models. A binary tree structure. 



obtained by hierarchically partitioning the original volume of size io in subvol- 
umes of size in = 2~"'£o, is used to describe fluctuations at different scales (fig. 
la illustrates the one-dimensional case). The energy dissipation, e„, associated 
to a cube at scale is multiplicatively linked to the energy dissipation, e„_i, 
at the larger scale, in~i, through a random variable Wn : 

e„ = Wnen-l = WnWn-lWn-2....Wie. (5) 
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The WnS are identically and independently distributed positive random vari- 
ables. The structure functions are now defined as 

S^f ^ E{el/^}Fj' (6) 

Using (^) and we obtain: 

= with Cip) = p/3 - log2 E{WP/'}, (7) 

where £'{■•■} denotes the mathematical expectation. Here, —\og2 E{Wp^^} 
is the correction to the K41 exponent in the structure function of order p. 
Actually, the multiplicative model has very interesting properties when 
considered in the light of the theory of large deviations, as noted, in particular 
by Oono and Collet and Kukiou It follows from Cramer's [|rB| (see 



also Refs. [|T^ and |2^) work that, for large n's, the quantity: 



- logs en = - (logs Wi + logs W2 + --- + logs Wn) (8) 
n n 

can deviate from its limit value, i?{logs W}, by a finite non- vanishing amount 
a with a probability which decreases as exp(— /(Q;)n). The CramerQ function 
/(a) is positive and convex and can be identified with an entropy in applications 
to statistical thermodynamics . The presence of such large deviations, with 
probability decreasing exponentially in n (i.e. as a power law in in), introduces 
fluctuations in the effective scaling exponent and therefore an overall change 
of the scaling properties of all structure functions. Indeed, we can rewrite (^ 
in the following form : 



As is well-known, in large-deviations theory, the function ({p) is given by a 
Legendre transformation : 

C(p) = inf (^(l-a) + /(a)/ln2). (10) 

Our aim now is to construct a deterministic variant of such models. We shall 
do it in two steps. First, we consider, instead of independent random variables 
{VFj}, successive points on the orbit of a Markov process. This means that we 
consider a Markov process on a phase space X with a transition probability 

name suee ested by Mandelbrot 
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operator P and an observable /?. : X — > R^. The (generalized) structure 
function at the scale in = 2^"'£o is then defined as : 

S^'^ = Slr\h,P] = El^f[^h^{x,)Y (11) 
where Xk & X are points of an orbit of the Markov process. Introducing now : 

n 

Z^^^Uh'ixk), (12) 

k=l 

we easily find that the conditional expectation of Zj^p) with respect to the 
initial distribution density, p(a;), is given by: 

E{Z^\P) = i hnx)p{x)Pl^ l{x) dx. (13) 

Here, the operator Ph^p is defined by the relation Ph^p(f){x) = h^{x)P(f){x). This 
concludes the first step. The second step is to change from a Markov random 
dependence to a deterministic dependence by use of a chaotic map. It seems 
that this can be done in an obvious way, replacing the transition probabil- 
ity operator by the transfer-matrix (Perron-Frobenius operator) of the chaotic 
map. However, some additional work is needed to define a deterministic con- 
struction possessing the same tree structure as in the random case (fig. la). 
Indeed, with a deterministic map, how can we avoid giving the same value 
to the two off-springs of the next generation, thereby trivializing the whole 
tree-structure? 

It is necessary to define a consistent procedure to distinguish the branches 
of the tree. We build up this branching-deterministic process by inserting a 
small amount of noise at each node and by considering the total process as 
the superposition of the deterministic transfer along consecutive levels plus the 
noise. The final map will be obtained by taking the zero- noise limit. 

Let f : X X E he the deterministic map which describes the relation 
between any two consecutive scales of the tree structure (see fig. lb), and let 
us fix an observable /i : X — > R^. Then for any fixed number r] > we may 

consider a random Markov chain (x^**^) such that 

4?i = /(4''^)+^e., (14) 

where the ^fc's are random variables, independently and identically distributed 
on the interval [—1, 1]. Now, we may apply the previous construction by consid- 
ering different realizations of the noise on each link connecting consecutive 
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nodes. Then, any cascade process along a branch of the tree in fig. lb will be 
described by the map (|T^) with different realizations of the noise. 

We define our structure functions as the zero-noise limit of (|TTp. 
For a fixed value of 77 > it follows from ([13|) that the large-n limit of the 
structure functions is governed by the spectrum of the operators : 

Ph,p,r, = h^Q^P. (15) 

Here, Qr? is the transition operator for the random perturbation, P is the 
transfer-matrix, or Perron-Frobenius operator of the map / and h is an ob- 



servable (for definitions see for example Ref. |p2|| ). 

Let us mention, incidentally, that there is also another possibility to define 
a deterministic (chaotic) process on a tree structure. Using a method of finite- 
state Markov approximations of chaotic maps, proposed by Ulam, we obtain 
a Markov chain with transition probabilities given by Pij = \Xi U f~^Xj\/\Xi\, 
where the set of Xj's defines a suitable finite partition of the original phase 
space X. It is known, for instance, that for sufficiently "good" maps (for 
example piecewise expanding with derivatives larger than two), invariant dis- 
tributions of this simple Markov chain converge to the density of the invariant 
measure of the map under consideration when |Xj| 0. With this Markov 
process, we construct the analog of (|T3|), viz the operator P^^p = h^P. (Here 
P is the transition operator with matrix elements Pij.) These possibilities will 
not be explored further here. 

We now observe that the construction based on taking the zero-noise limit 
appearing in (|T5|) is similar to that of the mathematical theory of fast dynamos 



23| , P^ , P5[] . Fast dynamo theory describes the phenomenon by which rapid 
magnetic field growth can be sustained in the presence of a prescribed velocity 
field, when taking the zero-diffusivity limit. From a formal point of view, fast 
dynamo theory involves a combination of two operators : a transfer-matrix 
for some deterministic map, associated to a deterministic velocity field, and a 
diffusion-like operator, or equivalently small-amplitude noise. 

The main purpose of the theory is then to find the properties of the zero- 
noise limit of the combined operator. Oseledets , in his study of the dynamo 



problem, considered a noise-perturbed Perron-Frobenius operator similar to 
that of ([l5|) and was able to show in special cases that it converges for r/ ^ 
to the Perron-Frobenius operator for the deterministic problem. The latter is 
than said to be "stochastically stable" . 

Returning to the structure function of the multiplicative model, we shall 
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also assume that the Perron- Frobenius operator is stochastically stable. In 
addition, we assume that the map (p!4D , with 77 = 0, is chaotic and ergodic. We 
can then compute the structure functions as averages along-the-orbit for the 
deterministic map, i.e. 

-r N n 

with ,Xj+i = f{xi). 



3 A class of deterministic cascade models based 
on shell models 

The construction defined in Sec |l|, after taking the zero- noise limit as explained, 
amounts simply to the following : we keep (|^) as it stands, but instead of 
searching the PVn's randomly, we assume that : 

W„ = ^7 W„_2, Wn-r) , (17) 

where g is a. deterministic map which involves a finite number r of antecedents. 

There is an alternative formulation in terms of characteristic velocities asso- 
ciated to the hierarchy of scales /„ = By analogy with @, we introduce 
a set of velocity variables, denoted «„, related to the e„'s by: 



\u. 



|3 



7 ^ ' 

The u„'s can be thought of as (real or complex) velocity amplitudes associated 
to eddy-motion on scale In- Instead of (H) we then use 

Un = qnqn-iqn-2----qiUQ, (19) 

while (^TD becomes : 

qn = fiqn-i, ,qn~r)- (20) 

Eq. (^) will be called the "ratio-map" since the g„'s are the ratios of successive 
velocity amplitudes. 

We shall now consider a particular class of ratio maps generated from "shell 
models" . The latter can be viewed as the poor man's Navier-Stokes equations : 
instead of the whole velocity field, one retains only a discrete set of velocity 
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amplitudes «„, the dynamics of which is governed by a set of coupled differential 
equations of the form : 

where the F„'s involve finitely many neighbors (in scale) of The various n's 
are referred to as "shells". Except for viscous and forcing terms, the F^s are 
chosen quadratic and satisfy energy conservation. We shall not here attempt 



any review of the considerable literature on shell models (see e.g. Refs. [|TOi |T2 



2B, E^, and references therein). Our intention is to use a static form of shell 



models, i.e. we assume : 

Fn{Un-s,---,Un,...,Un+s)=0- (22) 

It is not our intention to discuss such thorny issues as : are shell models "good" 
approximations to true turbulence? Can we learn something about the dynam- 
ics of shell models by studying the time-independent solutions? 

We just observe that static shell models immediately generate ratio-maps 
of the form (|20|). Indeed, ignoring viscous and forcing terms, which are not 
relevant in the inertial range, and assuming that Fn is a homogeneous polyno- 
mial of degree two in the m„'s, we find that g„ = Un/un-i satisfies a recursion 
relation of the form (pOf). An example will make this clear. One popular model 
is the Gledzer-Ohkitani-Yamada (GOY) model [|rT|, ^ which, in its complete 
form (with viscosity, forcing and time-dependence), reads: 

Ki + ^klu*^ + fn = Fn = -ikn (^Un+lUn+2 - ^Mn+l^n-l - ^Un-2Un-lj ■ (23) 

Here, z/ is the viscosity, /„ is the forcing, and fc„ = is the shell-wavenumber. 
The inviscid, unforced and static GOY model gives : 

Un+lUn+2 - ^Un+lUn~l - ^Un-2Un-l = 0. (24) 

4 o 



Hence, the recursion relation for successive ratios is 



1 



C[n-lC[n-2 



11 1 

4 8 qn-iqn-2qn-3 



(25) 



It turns out that the ratio-map generated by the GOY model is quite trivial. 
Indeed, from (^) it follows that the product Zn = Qn Qn-i 9n-2 satisfies a first 
order recurrence relation which has a single stable fixed point Zn = z^^ = 2~^. 
Hence, there is no chaos, and K41 scaling is obtained. To obtain chaotic 
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variations (in the shell index), we need to have more terms in One simple 
way is to perturb the GOY model by putting an admixture of another popular 



shell model, the Desnyanky-Novikov (DN) model |TI] 
the following ratio map : 



In this way we generate 



Qn — f{<ln-l,<ln-2,<ln-3), 

f{x,y,z) = ^((^ + l±-]+^(^-2y + 6{--2xy)]]. (26) 
xy \\4: SxyzJ \xy x J J 

Here, S is the so-called ratio of backward to forward cascade amplitudes in the 
DN model and 7 is the admixture parameter, which in subsequent numerical 
computations will typically be taken small (about 10~^). The hybrid ratio-map 
(p6D , which will be referred to as GOY-DN in the sequel, has still a K41 fixed 
point Qn = q* = 2"^/^, but it is easy to find a region in the phase-space of 
the parameters 7 and 6 where the K41 fixed point is unstable; for example 
7 > 0, V5. Further information about the behaviour of (p6|) requires numerical 
simulations. 



4 Numerical results 

We found numerically that the GOY-DN map, defined by (^) develops chaotic 
behavior for a sizable set of values of the control parameters 7 and 6. Most 
of the numerical simulations reported hereafter were done with 7 = 0.01 and 

5 = 4. However, the qualitative feature of the results appear to be very stable 
against changes in the parameters, as long as the dynamics are chaotics. 

Fig. 2a shows a typical orbit of the map (p6D. The initial condition is ob- 
tained by adding a very small random perturbation (r.m.s. value about 10~^) 
to the K41 fixed point value. The orbit looks chaotic and also displays strong 
intermittent fluctuations. Its largest Lyapunov exponent is found to be about 
0.22. The Kolmogorov fixed point, although it ceases to be stable as soon as 
7 > 0, still plays an important role in the dynamics. Fig. 2b shows an en- 
largement of fig 2a, revealing that the orbit mostly oscillates close to the K41 
fixed point, occasionally going on wild excursions. We mention that in the un- 
perturbed GOY map, the K41 fixed point is neutrally stable : two eigenvalues 
of the derivatives matrix, calculated at the fixed point, have modulus one. Of 
course, the fixed point becomes completely unstable as soon as 7 > 0. 

We now turn to the structure functions (|16]), evaluated by averaging along 
the orbit with several million points. Fig. 3 shows the structure functions for p 
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from one to four. They are seen to follow power-laws (in A;„) at large n's. Ex- 
cept for some residual chaotic noise, the exponents are exactly given by their 
K41 value: ({p) = p/3. Here, a word of warning is required. In numerical 
simulations of the full time-dependent shell models, it is not practical to have 
more than, say, 30 shells, since time steps decrease exponentially with n. Had 
we used such a low value of n for estimating exponents of structure functions, 
we would have predicted erroneous values, which would be misread as "mul- 
tifractal" corrections to K41. We do not want here to open the Pandora-box 
of whether the multifractal corrections detected in the time-dependent sim- 
ulations of shell models p6| , |27[] are or are not finite-shell artefacts (similar 



questions can be asked for finite- Reynolds number turbulence data). 



The presence of exact K41 scaling suggests that the GOY-DN model has 
no large deviations. This can be checked directly by plotting the p.d.f. of the 
normalized sums of logarithms : 

1 " 

S(n) = -^log2|gfc|, (27) 

playing in the deterministic case the role of (|^) in the random case. 

Fig. 4 shows this p.d.f. for two values of n. Note that as n grows it is 
increasingly peaked at the K41 value —1/3. It is thus likely that for n oo, 
it becomes a delta function at this value. A possible interpretation of this 
phenomenom is presented in the next Section. 



5 Neutrally unstable and sporadic maps 

We present a very simple example of a class of maps which displays chaos and 
no large deviations. Let us define : 

/^(x) = Frac {x^ + x) : [0, 1] [0, 1] , z> 2. (28) 

The deterministic maps are chaotic, piecewise expanding, and have a neu- 
tral (or indifferent) fixed point at the origin. These maps are "sporadic", in 
the following sense : The growth of infinitesimal errors is controlled by a 
stretched exponential with exponent less than one. Hence, the Lyapunov expo- 
nent is zero, although the maps certainly deserve to be called "chaotic". The 
derivative of any of these maps is equal to one at the origin and is strictly larger 
than one at all other points where it is well-defined. This leads to some unusual 
ergodic properties, for example the nonintegrability of the invariant measure. 
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the density of which has a power-law singularity near the origin [p8| , |29| , |30 
The reason is that once the orbit gets very close to the neutral fixed point, the 
time to escape from it is inversely proportional to a power of the distance from 
the fixed point. As a consequence most usual ergodic averages such as ( pT]) 
will be dominated by the contribution from the fixed point, so that there are 
no large deviations. 

As for the absence of large deviations in the GOY-DN map (^) we are led, 
at this moment, to speculate that something similar happens: the presence of 
a neutrally unstable fixed point (or, more likely, a neutrally unstable periodic 
orbit) dominates the large-n asymptotics of structure functions. 

Let us finally make two remarks. First, the positivity of the largest Lya- 
punov exponent in the GOY-DN map is not inconsistent, as one would 
think at first, with the kind of behavior displayed by the maps (^) which have 
zero Lyapunov exponent. To illustrate this, consider the direct product of two 
one-dimensional maps : the first is any map chosen in the class (^) and the 
second is the piecewise linear map g{y) = Frac(2|/) from [0, 1] to itself. The 
product has as its largest Lyapunov exponent log 2. Nevertheless, for the ob- 
servable h{x,y) = xy, ergodic averages are still dominated by the contribution 
from X = 0, as if the multiplier y was absent. Second, let us mention that 
for the maps (^) it is possible to give an estimate of the average fraction of 
time during which the dynamics is chaotic I^Oj. This average goes to zero as 
a power-law of the number of iterations n (e.g. as log n/n when z = 2). If 
the analogy between the sporadic maps and the GOY-DN model is correct, 
this suggests a possible interpretation of the rather strong corrections to K41 
scaling at finite n, which eventually disappear at large n. 
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Figure Captions 

Figure 1 : a) A one-dimensional representation of the dyadic tree for the ran- 
dom multiplicative chaos model. Each variable {Wi} is chosen independently 
and identically distributed, b) The branching process for the noisy determin- 
istic maps; the ^'s are chosen independently on each link. 

Figure 2 : a) A typical solution of the difference equation (pBI). b) A zoom of 
the previous figure. The horizontal solid line corresponds to the Kolmogorov 
fixed point. 

Figure 3 : Log-plot of the structure functions S^p^ as a functions of n, for 
p = 1,2, 3, 4. The solid lines correspond to the Kolmogorov 1941 scaling. 

Figure 4: Probability density functions (p.d.f.) of S(r;,), defined by ( pT] ) 
for n = 21 and n = 36. Notice that the peak at the Kolmogorov value, 
E(ri) = — 1/3 becomes more pronounced as n increases. 
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